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1 Introduction 



Light cone quantization was introduced by Dirac in his exploration pP of possible repre- 
sentations of the Poincare Group for relativistic formulations of quantum mechanics. The 
most obvious application in nuclear and particle physics is to high momentum transfer pro- 
cesses for hadrons and nuclei, since in a light cone representation a diagonal Lorenz boost 
operator can be defined, while in the standard instant form representation Lorentz boosts 
of bound systems are difficult to formulate. 

There have been several reviews of light cone quantization for Hamiltonian Dynamics 
with various applications to bound states (see Ref. [2] for references) and the light cone 
representations of bound states in supersymmetric theories|2]. The light cone community 
holds an annual meeting, and this review is based in part on developments in this area 
presented and discussed at the International Workshop Light cone 2002 j4j 

The light cone formulation of Quantum Chromodynamics (QCD) field theory, the presently 
accepted theory of strong interactions, is particularly important, since only by treating pro- 
cesses from low momentum transfer to high momentum transfer can one explore the transi- 
tion from nonperturbative to perturbative regions, which is still an unsolved problem. Light 
cone formulations are important for such studies. E.g., in one of the early applications of light 
cone formulations to form factors it was shown that with a model that predicts that non- 
perturbative processes dominate the pion form factor until momentum transfer ~ AGeV^, 
and that at > 5 GeV^ perturbative QCD theory might begin to be applicable, it was 
found in an instant form treatment of the same model [HI that only at much higher values 
of can perturbative QCD be used. Light cone treatments of the pion form factor are 
reviewed below. 

In the present review we shall center on applications to processes involving either high 
momentum transfer or both high and low momentum transfer. 

2 Light Cone Quantum Mechanics 

In this section we review the light cone formulation of relativistic quantum mechanics. 
Quantum mechanics is based on states and operators, which must tranform correctly under 
Lorentz transformations in a relativistic theory. The Poincare group of translations, rota- 
tions, and Lorentz transformations in the standard instant form has serious problems for 
Lorentz transformations, and as we shall see this problem is much less serious in the light 
cone representation. 

These problems are more serious in field theories, which we review in the following section, 
and a number of the concepts which arise in quantum mechanics apply with modifications 
in field theories. 
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2.1 Lorentz Transformations and the Poincare Group 

The invariant distance, ds^ is defined in terms of the four position, by the metric tensor, 
by 

ds" = g^^x^x", (1) 

where the Greek indices run (1,...4). An inhomogeneous Lorentz transformation of the 
position 4- vector is given by 

x'^ = a^ + A'^x", (2) 
where a^^ is a constant and is the constant Lorentz matrix, which satisfies 

A^A^^ = (3) 

The ten generators of infinitesimal Lorentz transformations are: 

P'* are four momentum operators (4) 
M^'^ give six independent angular momentum operators 

with Mf"" = -M^t". 

These ten operators form the Poincare group, which satisfy the the commutation rules 
(Poincare Algebra) 

[P'^,P"] = (5) 

[M^'^M^"] = iig^'^M"^ + g''^M^"' - g^'^M"" - g'^'M^'^) 
In the next two subsections we consider the instant form and light cone representations. 

2.1.1 Instant Form for Dirac Particles 

The instant form is the one commonly used in relativistic quantum mechanics: 

= (x",x) = {t,x\x^,x^ ^ z) (6) 
g^'^ = (1,-1,-1,-1) 

g^"" = Ofor^T^i/, (7) 

with fjL = (0,1,2,3). 

Let us first consider a free Dirac Particle. The ten generators (c= h/27r =1) at t=0 are 

P" = (P°,P) (8) 
Ji = eijkM^^ = {fxP)i + Si 
Ki — Moi — (7oi, 



where cr^i, = ^[7^, 7j,]/2. and S are the usual momentum and internal spin operators. 
Wave functions transform under the group by 

translations ip'ix + X) = e'^'^ipix) (9) 
rotations ?/^'(f^ = e^^^^->(x) (10) 

Lorentz transformations ip'{x^') = e'''^^"^ip{x^), (11) 

with a rotation of 9 about the n axis and a Lorentz boost of tanhcu = v in the n direction. 

For a system of interacting fermions in the instant form the boost operator K contains 
the interaction potential. There has been a great deal of literature on this subject, starting 
with Foldy's expansion in c~^ We do not discuss this further here. 



2.1.2 Light Cone Form for Dirac Particles 

The light cone form is defined by 
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In this light cone form one can choose seven generators to be independent of the inter- 
action (often called "good" generators), such as (see [HElin]) 



Good: P+, P\ P^, El, E2, Ks, J3, 



(16) 
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so that one can carry out a Lorentz boost in one direction without involving interactions. 
This is also true for field theory, as we shall see, and is most significant, since interactions 
generally produce particles. The other three generators (the "bad" generators) 

Bad : Fi, F2, P" (17) 

contain the interactions. Therefore, although a boost in the z direction is interaction free, 
rotations about the x and y axes are not; therefore prescriptions for total angular momentum 
for composite states are needed. See, e.g., Ref[TU]. 



2.2 Light Front Hamiltonian Dynamics 

One method for constructing a relativistic quantum theory for interacting particles with an 
interaction-independent Lorentz boost follows the Bakamjian- Thomas idea^T] of using M, 
the invariant mass operator. See Ref.fT^ for a review and references. 
The energy operator in this method is P~ 

P - (IB) 

The method of obtaining solutions to the Hamiltonian eigenvalue problem, with the Hamil- 
tonian operator H containing the interactions, 

> = — 1^> (19) 

for the states and wave functions by a Fock state expansion is discussed in detail in Ref.j2], 
with discussion of applications to QCD and QED. There is an extensive discussion of ap- 
plications for relativistic corrections in nuclear physics in Ref.|12j. In the present review we 
concentrate on applications of hght cone field theory, which is discussed next. 



3 Light Cone Field Theory 

In this section we review the light cone formalism for quantum field theory. One of our main 
objectives is to discuss the derivation of the ten Poincare generators, with an interaction-free 
Lorentz boost operator. First we review the instant form procedure for deriving the Poincare 
generators. 

3.1 Poincare Generators in Instant Form Field Theory 

The starting point in a quantum field theory is the Lagrangian density, C{x) = C{tp'^, d^ip°'), 
with the fields ip", the conjugate momenta 7r^(a;, t) = dC{x)/{dQ'ipf3{x,t)), and the commu- 
tation rules of the fields, 

[V^"(f,t),7r^(f,t)] = t5^,f35ix-x'). (20) 
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The Poincare generators are obtained from the energy momentum tensor, which can be 
derived from 

E.g., the Lagrangian density for pure glue QCD is 

_ \g-G (22) 
4 

with 

G^, = d^A,-d,A^-ig[A,,AJ\ (23) 
A, = AIX^I2 

where A" are the eight SU(3) Gell-Mann matrices, ([Aa, A;,] = 2ifabc^c)- The corresponding 
energy momentum tensor is given by 

Ti9iue),u ^ Y^^G>:,-G:,,-\g^^GfG^pa). (24) 

a 

The momentum operators are 

pt'{x) = T^i'ix) (25) 

and the Hamiltonian density is 

n{x) = T°°. (26) 
The form of the M^^ tensor depends on the specific theory. We consider two examples next. 

3.1.1 Poincare Operators for Scalar Fields 

For a scalar field with the Lagrangian density 

C ^ l(a^0(x))2-l</,(x)2 + Ant(x), (27) 

one finds 

W = T" = (a„4'(x)f - £i„,(i) (28) 



where the boost operator has been evaluated at t=0. It is clear that the contain the 
interactions. 
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3.1.2 Poincare Operators for Dirac Fields 

The Larangian for interacting Dirac particles is 

C{x) = |[V;(x)7'^9^V'(^)-(9^V5(^))7'^V'(^)]-W(x)V(x)+Ant(^)- (29) 

One finds that the internal spin leads to an additional term in the angular momentum 
operators: 

/?■ f)C 
d^xhc^'T^i' - xi'T^'' - , (T^>] . (30) 

2 d{dQil)) 

Using — M°', at t=0 one obtains for the boost operator 

K ^ j d'x[]^xn{x) - |v(a^)7V'(^)] ■ (31) 

As expected, the boost operator contains the interactions. Thus the state obtained by a 
boost with velocity v = vn the state I'V'o >; 

|V.,> = e^-^-"|V'o>, (32) 

will be quite different form the state at rest, since the boost K contains fields which create 
and destroy particles. This is a serious problem for studies of, say, form factors of hadrons 
at high momentum transfer. 



3.2 Poincare Generators in Light Cone Field Theory 

Since the Poincare generators depend on the specific Lagrangian, as discussed above, let 
us consider the special case of the light cone representation of the Poincare operators for 
the scalar field theory, .x^ = (a;"*", x\ x^, a;~) = (a;"*" , a;_L , x~ ) and C = |(c)^0(a;))^ — |0(a;)^ + 
'^int{x). The field satisfies the commutation rule [(f){x), 4>{y)] = S{x± — y±)S{x~^ — y~^)S{x~ — 

The energy momentum tensor has "good" operators 

T++ = (9+0)2 (33) 
= d+(j)d,(j) ^ = (l,2) 



and the "bad" operator 
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From this we obtain the seven "good" Poincare generators (i=(l,2)) 

P+ = J d^x^dx-{d+(t){x)f (35) 

= J d^x^dx~d+(j){x)di(j){x) 

= J d^x±dx-[x\d+(p{x))di(j){x) - x'^{d+(f){x))d2(f){x)] 

= j d^x^dx-[x-^{d+(p{x)f - x^{d+<p{x))d2(t){x)] 

M+' = I d^x^dx~[x+{d+(t){x))di(t){x)-x\d+(t){x)f]. 



Again we find that the boost in the z direction is interaction-free. The three "bad" Poincare 
generators, containing the interaction, are 



P- = 



H = j d^x^dx-[^iV±(l)f + m^(j)^)-Cint\ (36) 

= j d^x^dx-[x-{d+^{x))d4{x)-x\]^{y^(t)f + m^'')-C,nt)]- 

This completes our review of hght cone field theory. We now review the light cone 
formulation of the Bethe-Salpeter (B-S) equation, from which one obtains the B-S ampli- 
tudes needed for all hadronic studies, and the Schwinger-Dyson (S-D) equation, from which 
one obtains the dressed quark propagator needed for the B-S equations for all mesons and 
baryons. 

4 Bethe-Salpeter in a Light Cone Representation and 
the Pion Form Factor 

In this section the Bethe-Salpeter (B-S), Schwinger-Dyson(S-D) formalism is reviewed. Only 
quark-antiquark systems are considered for the B-S equation, with application to the pion 
form factor. We refer to a quark-antiquark system as a two-quark system for simplicity. In 
the present section we review applications to the pion form factor. In the next section a 
recent light cone solution of the quark S-D equation is discussed 

4.1 Instant Form of the Bethe-Salpeter Equation 

The Bethe-Salpeter equation is an exact equation for the two-particle propagator. For a 
derivation see, e.g., ReffT^. For a two-quark propagator the B-S equation is 



m - S)(i - m - S)S'(xi,X2;?/i,?/2)2<7 = S{xi,yi)S{x2,y2) (37) 



+ 



j d'^Zid'^Z2K{xiX2] ziZ2)S{zi, Z2\ 2/1, 2/2)29, 
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Figure 1: Schematic form of B-S amplitude 



where is the two-quark propagator, S is the quark propagator, m is the current quark 
mass, S is the quark self-energy (from the S-D equation, shown below), and K is the B-S 
kernel. 

If the two-particle system has a bound state with mass M, then there is a pole in S2q in 
the total momentum P variable, 

S{xi,X2;yi,y2>2q\p2^M2 = ^ p2 ^2 • y^^) 

The quantity \E', which in the rest system of the bound state of the two-quark system \M > 
is 

^(a;i,X2) = <0|T[^(a;i)^(x2)]|M>, (39) 

which is often referred to as the B-S wave function or B-S amplitude. It satisfies the B-S 
equation 

{i -m-T){i - m - S)^(xi,X2) = j d'^yid^y2K{xiX2;yiy2)^iyiy2)- (40) 

The B-S amplitude and equation is pictured in Fig ^ Note that the solution requires 
a set of coupled equations for the vertices, F, self mass of the quark (S-D equation), and 
the B-S amplitude. The Kernal for a perturbative calculation is given by the sum of all 
irreducible diagrams. For nonperturbative theories like QCD one must use models. 



4.1.1 Problems With Instant Form of B-S Equation 

About a half century ago an important paper by Wick^] pointed out a number of serious 
problems with the B-S equation and B-S amplitude. One of the problems is the dependence 
on the relative time (or relative energy) variable, whose physical meaning never was clear. 
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He found a solution to this problem via the Wick rotation to Euclidean space. This has 
wide application, and is still used in hadronic physics. Since this problem has been solved 
by the light cone formulation we do not review this idea, but aspects will appear later in 
this review. 

4.2 Light Cone Bethe-Salpeter Equation 

As discussed above, in order to obtain a B-S wave function for hadronic physics one must 
also solve the Schwinger-Dyson equation for the dressed quark propagator, which we discuss 
in the following section. In any case, models are needed for QCD. Here we discuss some 
early models of the light cone formulation of the B-S equation and applications to the pion 
form factor. 

Since in the infinite momentum frame the pair creation and annihilation processes are 
suppressed |2Uj. it has many of the attributes of the light cone formulation. In an early 
paper [21] QED at infinite momentum was treated, and this work has been used recently for 
obtaining a light cone S-D equation, discussed in the next section. The light cone form of the 
B-S equation can be obtained by using the cxD-momentum frame or projection onto the 
light cone |23j. In momentum space, using the notation for the momenta of the two quarks 
ki = (xiP"*", k±, ki), ki = {x2P~^, —k±, fc^), with k~ = {w? + k\)/kf , the light cone equation 
for the B-S wave function, V'l^i; ^2, P) is 



In this form one assumes that the mass, m, of the quark is a constant constituent quark 
mass, and uses a model for the kernel. 

4.3 Application of Light Cone B-S AmpHtudes to Pion Form Fac- 



In this subsection we consider applications of the light cone formalism to the pion form factor. 
The importance of the light cone formalism to study the transition of the pion form factor 
from a low-momentum region in which nonperturbative effects dominate is first reviewed, 
followed by treatment with a fixed quark mass, and then the most recent work with a running 
quark mass. See Sections 8 and 9 for details on 1-c wave functions. 

4.3.1 Light Cone Wave Function and Transition to Asymptotic QCD 

One of the interesting problems in the past two decades is the nature of the transition from 
the low-momentum region of nonperturbative QCD to the asymptotic region of perturbative 
QCD. This problem has been studied theoretically for many years, such as the modified 




(41) 




tor 
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quark distribution functions j21], and experimentally. An early application of light cone 
methods 123 was simply the use of a light cone vs instant form pionic wave functions. 
The result is quite startling and a dramatic illustration of the advantage of a light cone 
representation. 

Given the wave function of the pion, 0^, for the low-momentum or soft part of the pion 
form factor with a nonrelativistic wave function the pion form factor is simply 

Fr^'iQ') = Jd'p<j>:{p,P + Q)Mp,P). (42) 

On the other hand it is known that the high momentum transfer (hard) limit perturbation 
diagrams can be used, with the result j^H] 

FAQ')\q^^oo ^ Snf^^^, (43) 

with a similar result obtained by others^22|. It is of great interest to know when perturbation 
theory can be applied, which is when the nonperturbative soft contributions become smaller 
than the high-Q contributions. In an attempt to study this a nonrelativistic harmonic 
oscillator wave function 

MP) = jj^.^''"^^"'-' (44) 

was used|6|. Even at ~ 10 GeV^ the soft part dominates the hard term. On the other 
hand, recognizing that even at 1 GeV^ relativistic effects are important, a calculation with 
the the same wave function in the light front representation jj]: 

0^(p,P) = A^e (45) 

with N a normalization constant, has been carried out. The results are shown in Fig |21 
As one can see by comparing the light cone (LC) with the instant form (IF) curves, at 
greater than about 0.5 GeV^ the relativistic corrections become increasingly important. The 
curve labeled ^ oo gives the limit of the hard scattering processjSEI- light cone 

model the soft part of the form factor becomes smaller than the asymptotic hard part for 
Q"^ ~ AGeV"^, while in the nonrelativistic instant formP^ the soft part dominates until very 
high Q^. A vastly different interpretation of the transition from soft to hard QCD results. 

4.3.2 Light Cone B-S Amplitude for Pion Form Factor 

There have been a number of studies of the soft (nonperturbative) pion elastic form factor 
using hght cone or other relativistic models j2Hl 1211 and an early studyjHO] using light cone 
wave functions to test the amplitudes of Ref . [21] , as well as the work of Ref . . 

In this subsection we discuss the study of the B-S equation for the pion amplitude [HT] 
with a calculation of the pion form factor including both the soft and hard parts. The 
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Figure 2: Pion form factor in a light cone representation (LC, solid) and with an instant 
form (IF, dotted) harmonic oscillator wave function 



starting point is the light cone B-S equation, Eq. ()4Hl . using the kernel K = Ki + K2, with 
K2 the one gluon exchange kernel used to get the hard asymptotic form for the pion form 
factor |2S1 EZI, and Ki modeled to give the soft confining part: 



K{xi,X2,k±,yi,y2,h 



Xl. 



xiX2{yi - y2y 



Xl 



y2 



(46) 



+ (l--2). 



(2/1 - xiY' 

with b an interaction strength, and m is a mass parameter. The tt"*" form factor is given by 



dx d?k±,2 
1 - x2 167r3 ^3 



-q,x) + U>\k^~^ 



X 



g,x)]^(A?j 



(47) 



The B-S amplitude is constrained by 1) normalization, 2) the pion decay constant, 
and 3) the pion mass. The pion mass is quite difficult to fit with this model and the 
best fit to the other parameters resulted in a pion mass of 356 MeV. To obtain a solu- 
tion, the spinors are first projected out. Expressing the B-S amplitude as \l'(xi, 0:2, /c_l) = 
$(a;i, X2, k±)V{pi)V{p2)XF, where V{p) is the melosh operator [221 CHI and Xp are the quark 
light cone spinors. After projecting out the spinors, the B-S equation for $(xi,X2,A;x) is 
solved using the technique of expansion in hyperspherical harmonics |22j. 

The result for the pion form factor is shown in Fig. 01 The prediction is that the 
asymptotic part of the B-S amplitude dominates by 10-15 GeV^. Other light cone model 
calculations [23 gave similar results. 
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Figure 3: Pion form factor in from B-S amplitude calculated with a light cone kernal con- 
taining confining and asymptotic parts 



4.3.3 The Running Quark Mass and the Pion Form Factor 

In addition to the transition from soft to hard processes dominating hadronic form factors, 
as one goes from the nonperturbative to the perturbative regions of momentum transfer the 
effective quark mass also undergoes an evolution from the constituent mass of low-energy 
quark models to the current quark mass in the QCD Lagrangian. This can be understood 
most directly by considering the quark propagator, which has the form 

Sip)-' = Aip') ^-Bip'), (48) 

where the running quark mass is given by M{p^) = B{p^)/A{p^). The functions A{p^), B{p^) 
are determined by a set of self-consistent equations in the Schwinger Dyson formalism, which 
we review in the following section. M(p^) is expected to undergo its transition from the 
constituent mass of about 300 MeV to a few Mev for u and d quarks at about 1 GeV. 

The effect of this running quark mass on the pion form factor was studied some years 
agojSSl, and in greater detail recently jE]. The pion form factor, given by Eq. ljTTj) and the 
discussion following that equation, is modified if one includes the dressed quark vertex. The 
elastic pion form factor can be written as 



X ^r2i#*,.,,(..kj. (49) 



'Pa 



where P^=]5g^=(l — x)P^ and k'x=k_L — xq_L in the initial pion rest frame, P_l=0. The 
helicity of the quark(antiquark) is denoted as Ag(q). The B-S amplitude is the one used in 
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Ref.[SI] 



dk 

^<l>(a;,kx)7^A,A,-(x,kJ 



(50) 



where defined in Eq. 

function. The Ball-Chiu 

section) 



BC 



is the radial wave function and y is the spinor wave 
ansatz is used for the vertex (see the discussion in the next 
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-{p+p'r 



-r 



{p+p'r 



(51) 



p^ _ p^ 

In order to consistently determine the running quark mass, the S-D equation must be 
solved in a light cone representation, which is reviewed in the next section. In Ref. jjl] models 
were used. For the crossing antisymmetric (CA) parameterization the model is 

l + exp(-/iVA4) 



M{p ) = mo + {nic — mo) 



l + exp[(p4-;,4)/^ 



41' 



(52) 



where mo = 5 MeV and rric = 220 MeV are the current and constituent quark masses, 
respectively. The parameters /x and A are used to adjust the shape of the mass evolution. 
For the CA picture the two sets of parameters used are (/i^,A^) = (0.95,0.63) [Set 1] and 
(0.28,0.55) [Set 2] GeV^. The resuhing M{Q^), with = -p"^, is shown in Fig. H 

The results of this model for the pion form factor are shown in Fig. El In comparison 
with Fig. El it is clear that the running quark mass must be taken into consideration in order 
to trace the evolution from the soft perturbative region to the asymptotic region. 



4.3.4 Other Light Cone Calculations of the Pion Form Factor 

There has not been much experimental progress in extending accurate measurements of the 
pion elastic form factor to higher momentum transfers, and there is little hope in that the 
region above 10 GeV^ will be reached soon. There have been suggestions for improving the 
perturbative calculations, such as Refs. |HS1 EHl EZj , but detailed calculations of the pion form 
factor from the soft to the hard regions of momentum transfer require models of the nonper- 
turbative part, and little progress has been made over the past decade. This is unfortunate, 
since the soft to the hard transition for the pion occurs at a much lower momentum transfer 
than for the nucleon, but the experiment is quite difficult. 



4.4 Low Energy Theorems, Chiral Symmetry, and Light Cone 
Field Theory 

Although light cone field theory is generally most important at high momentum transfers, 
due to the difficulty in carrying out Lorentz transformations in the instant form (as discussed 
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Figure 4: The Running Quark Mass. Set 1, Set 2 defined in text 
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Figure 5: Pion Form Factor From a soft B-S Amplitude With An Evolving Quark Mass 



above), there has been at least one important application at very low energies. We briefly 
review the soft pion theorems for threshold pion photoproduction and the importance of a 
light cone formulation. A complete discussion with references can be found in Ref. ^j. 

The amplitude for pion photoproduction with a nucleon target, A^(7,7r)A^ has three 
isospin amplitudes, A^^\ A^~\ A^'^\ A model-independent expansion in the low energy pa- 
rameter /i = mT^/M^ gives 

A(-) = l + 0(/i2) 

= -^/2 + 0(/.2), (53) 

or in the soft pion limit 

A(-) = 1. 

= (54) 

Standard quark models have a great difficulty with these model-independent predictions. 
The main problem is that the pion photoproduction processes are related to pion exchange 
currents, and standard quark models give very different predictions at low energy than 
hadronic models jl^. As discussed in Ref. [21], the low-energy limit in standard quark 
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models gives 

= (55) 
9 

This is a serious violation of low energy theorems, which is also found in chiral quark 
modelsjlO]. 

This problem is resolved in light cone field theory models. An interesting observation 
is that for the calculation of A^~^ the standard pion exchangs process is replaced by the 
instantaneous term|27j. It was shown in jHEj that with a standard light cone harmonic 
oscillator model 

(j){x,k±) ~ eV ~^^) (56) 

the low-energy theorem predictions of Eq. ()53|) are satisfied. It should also be noted that 
the calulation of the Eo+ amplitude does not agree with low-energy theorems, and that 
the measurement of the p(7, 7r°)p threshold cross section [41^ obtained a result for the Ea+ 
amplitude also in disagreement with low-energy theorems. 

5 Light Cone Representation of the Quark Schwinger- 
Dyson Equation 

In this section we discuss the light cone formulation the quark Schwinger-Dyson Equation 
(SDK). There has been a great deal of interest in the SDE for hadronic physics, with excellent 
reviews 1^ . The Schwinger-Dyson (S-D) formalism is a technique for obtaining the 
dressed quark propagator; although this propagator is not in itself physical, it can be used 
to obtain nonperturbative hadronic strucdture. We review this in the next section. The 
quark propagator in space-time is defined as the correlator 

Sf. = <0|r[g;(x)gt(0)]|0>, (57) 

where T is the time-ordering operator. Since the propagator is diagonal in color, the prop- 
agator can be written as S{x) =< 0|T[g'^(x)g"(0)]|0 >, with no sum on the color index a. 
The S-D formalism starts with an exact expression for the quark propagator and explores 
models to obtain solutions in the presence of nonperturbative effects. 

5.1 Schwinger-Dyson Formalism 

The full dressed quark propagator satisfies the Schwinger-Dyson equation (SDE). (See FiglH)): 

s-\p) = s,\p)-i:{S{p)) 

S(5(p)) = j ^,l'^S{q)Tl{q,p)Df^{p-q) (58) 
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where r^(g,p) and D'j^^ijp — q) are the dressed quark-gluon vertex and dressed gluon prop- 
agator, Xa is a SU(3) operator in color space, 7^ is a Dirac operator, with Greek letters 
representing Lorentz indices and latin letters standing for color indices. The SDE is illus- 
trated in Figini 

The renormalized vertex is obtained from the vertex SDE 

r^(p,p') = Z7^ + |^y"^^'(p' + gK(p + g,p' + g)^'(p + g)i^V,P,9), (59) 

where Z is a renormalization constant and the kernel K{p',p,q) is related to the dressed 
gluon propagator D{q). This is illustrated in Fig. [7| 

The solutions of the SDE are of the form given in Eq. ()48p 

Sip)-' = A{p')^-B{p'). 

with the effective mass given by Mip^) = B{p'^)/A{p'^). 

The necessary ingredients for obtaining solutions are the dressed vertex and kernel. In 
practice, the dressed vertex is modeled using symmetries, and the main result of S-D cal- 
culations for applications to hadronic physics, in addition to the dressed quark propagator 
itself, is the dressed Gluon kernel. The transverse part of the dressed vertex is constrained 
by the Ward-Takahashi identity (see Ref|ISj for detailed discussions) 

{p — p)^T''{p,p') = S^'{p') — S^'{p) Ward — Takahashi Identity, (60) 

which is an extension to nonvanishing momentum transfer of the Ward identity, 

d 

T''(p,p') = — S(p) Ward Identity. (61) 

^P^^ 

This leads to the Ball-Chiu ansatz, Eq. (|KT|l . and the Curtis- Pennington forms jilj. discussed 
in detail in Ref.f^. 

The general form of the gluon propagator is 

where ^ is the gauge parameter. The Landau gauge corresponds to the choice of ^ = 0, while 
a so-called Feynman-like gauge corresponds to the choice 

D;1{x) = 6abS,M^) (63) 

for the model gluon propagator. 
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Figure 6: Diagrammatic representation of the SDE for the quark propagator. 
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Figure 7: Diagrammatic representation of the SDE for the dressed gluon-quark vertex. 
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5.1.1 Constraints on Quark Propagator 



One approach to the treatment of the nonperturative aspects of QCD correlators is to use the 
operator product expansion, which implies a continuation to large momenta. The operator 
product expansion of the quark propagator is 

<0|r[g;(x)g;:(0)]|0> = + + + (64) 

where < gg > is the quark condensate and < qa ■ Gq >=< qcy^vG^'^q > is the mixed 
condensate. This is illustrated in Fig. |S1 The condensates are vacuum expectation values 




Figure 8: Diagrammatic representation of the operator product expansion for the quark 
propagator. 



of local operators, and therfore are gauge invariant constants. The quark condensate has 
been known for many decades, and the mixed constant is somewhat uncertain, but has been 
estimated using the QCD sum rule method. These condensates can be used to constrain the 
S-D solutions by the relationships in Euclidean space 



< 0| : g(0)g(0) : |0 > 
< 0| : g(0)^a-G(0)g(0) : |0 > 



47r2 



+00 



dSS 



B{S) 



47r2 



+00 



dSS{ 



SA^{S) + B^{S) 
81B{S)[2SA{S){A{S) - 1) + B\S)] 



(65) 



16(5^2(5) + 52(5)) 



^ B{S){2-A{S)) 
SA^{S) + B^{S)^ 



with the notation that S = Q'^ the Euclidean squared momentum. Another constraint is 
given by 



•I IT 



47r2 



dS 



SM{S) 



A{S)[S + AP{S)Y 



M{S) 



SdM' 
2~dS 



(66) 



with M=B/A. 
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5.1.2 Rainbow Approximation For SDE 



The rainbow approximation corresponds to tlie clioice of 

= l^^al'^ rainbow approximation). (67) 



Tliis is the main feature of the Global color model j45], reviewed in Ref . j43j . Using the 
Feynman-like gluon propagator of Eq. ()63|) the self mass is given as 

^^^^ = ^^?s/^7'^^(g)7.^(g,p)^(p-g)- (68) 

In this rainbow approximation the quark SDE becomes a set of coupled integral equations 

d^q A{q^)q-p 



[A{p')-l] = 9llj^.D{p-q) 



q^A\q^)+B{q 



2^ 



There have been many applications of this model. The vertex for the SDE, Eq. ()59|) without 
the perturbative term has the form of a B-S amplitude for a bound quark-antiquark state, 
i.e., a meson. The B-S equation for a meson as a bound state with momentum P is 

rt{p,P) = 1 0^,S'^{q + xP)Tt{q,P)S\q-il-x)P)K{p,q,P), (70) 

where a,b are the color indices of the quark, antiquark, x is the momentum fraction of the 
quark, and K is the kernel. If one is given solutions to the quark SDE this provides a 
model of mesonic properties (see, e.g., [IHlllZlllHllini)- As discussed in the previous section, 
one of the most important studies of hadronic properties is the transition from the soft 
nonperturbative to the hard perturbative region of momentum transfer, which is best done 
in a light cone representation. For this reason it would be most valuable to have a light cone 
S-D solution to use to get hadronic B-S amplitudes, analogous to Eq. (|70|l . We now discuss 
the first successful light cone solutions of the quark SDE. 



5.2 SDE and Quark Propagator in a Light Cone Representation 

The method used for solving Eq. (j58|) in a light cone representation in Ref.pS) is to use 
a technique introduced originally for perturbation theory [HF- Chang and Ma showed how 
the rules of light cone perturbation theory [LCPT] can be derived in the infinite momentum 
frame from the usual covariant Feynman rules simply by changing into light cone variables: 

g±=g°±g3, q^={q\q^) 

An important feature of the new rules is that the range of integration over the q~^ variable 
becomes finite 

r"+00 fp'^ 

dq^ — > / dq'^, 
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where is the "plus" (longitudinal) component of the external momentum. This feature 
is related to the fact that in LCPT there are no diagrams with lines going backwards in 
time and no vacuum diagrams (with the exception of zero-modes and instantaneous terms in 
fermion propagators). This is a crucial property of light cone field theory, since it simplifies 
tremendously the structure of the vacuum. 

An important consideration is the form of the integrals needed for the calculation of the 
quark self-energy S in the SDE. We are faced with the non-perturbative self-energy diagram 
seen in Fig. The integrals occuring in the SDE diagram can be written as a sum of terms 
of the form: 

/ wy"^'^"^'' ■ ■ ■ ^''"^'^ ^ '''' " '^^^ ^ ''^ ' ^^^^ 

where the fq contains factors coming from the quark propagator and possibly from the 
dressed vertex function, fc contains factors coming from the gluon propagator and also 
possibly from the dressed vertex function, and the factors of g^' come from Dirac operators. 
Since the external momentum p is held fixed during this integration, possible factors of p'^^ 
and possible dependence of the Green's functions on are not relevant to the analysis below. 
Consider first the scalar case, = in ()71|) . i.e., no factors q^. Using the variables: 

a = q+/p+, s' = q-q, s = p ■ p, q'^ = q± - a p±, 

the integral in ()7H) becomes: 

/ + ''^ [-^'"n<iL s, a, e)) , (72) 

where all the integrals are over the entire real line, and 

V{ql, s, s', a, e) = ql + a{l — a){—s) + (1 — a)s' — iae 

In Ref . it was shown that if all the singularities of the functions fq and fa occur on 
or below the real axis, only the interval (0, 1) in the integration over a contributes to the 
integral. The model must ensure that the integrals are ultravioletly convergent and thus 
closing the contour with a semicircle at infinity introduces no additional contribution. The 
integral in (f72|) then becomes 

1 /'+°° ds' f (i'^~^q' / / \ 

4^ 2^/, "^"y 12^0^" '^'^ ^ '"^ [-a-'V{ql, s, s' , a, e) j . (73) 

It should be noted that to preserve important features of light cone theory (here reflected in 
the integral over a being over a finite range) one must make assumptions about the location 
of the singularities in fq and fc- These assumptions are much like the ones necessary to 
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justify a Wick rotation, however, the method of Ref.fTH] allows one to solve the quark SDE 
for time-like values of the momentum p, which is difficult to do after a Wick rotation. 
The model gluon propagator is 



-1 



k'^ + ie 



(74) 



The parameter x = for a Feynman-like gauge see Eq. (|H!?jl and x = 1 for Landau gauge. 
It has been shown [Slj that renormalization group arguments yield an approximate relation 
between the renormalized coupling constant, the renormalized gluon propagator and the 
effective couphng: 



glDnik' 



(75) 



where the subscript R denotes renormalized quantities, and Qeff is the effective running 
coupling constant. One model used the polynomial form 



N 



a 



eff 



1=1 



So 



k"^ + ie 



(76) 



and the parameters of the model were chosen to fit the PDG 52j values, shown in Fig. |H1 A 
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Figure 9: Fits to Oe// with model of Ref.jTSI a) for energies up to 350 GeV and b) a close-up 
for lower energies. 

sample result of the calculation using the iterative procedure is shown in Fig. 

In Fig. ^2 the results of a calculation of the light cone SDE are shown, with Se = Q^, 
the Euclidean momentum squared. The parameters for this calculation gave a reasonably 
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Figure 10: Comparison of the numerical and analytical calculations of the MS-renormalized 
quark propagator after one iteration in an ultraviolet divergent case, a) analytical subtrac- 
tion in numerical integrations b) with tail contribution and MS subtraction, c) without tail 
contribution or MS subtraction. 



good fit to and <: qq :>, but a rather large value for the mixed condensate. It has 
a strong infra-red enhancement. The results for the effective quark mass are too large at 
low momentum, but the evolution in is reasonable. A most interesting aspect of this 
calculation is the ability to calculate the quark propagator in the time-like region. 

The polynomial model used to obtain the results in Figs. IHl ^2 describes an infrared- 
enhanced gluon propagator. Although this can give confinement, it cannot give the QCD 
nonperturbative midrange attraction. On the other hand the instanton model can give the 
strong midrange attraction, even though it does not give the observed confinement. The 
instanton form of the gauge color field is obtained in a SU(2) Euclidean classical model 



192p4 

+ p2)4 



where rja^v are given in Ref.^j and for p, the instanton size, the instanton liquid model [55] 
is used. The quark propagator in the instanton/anti- instanton medium, S"/, has been derived 



24 



in this model jnSl; giving 



Slip) 

MP) 

BiiP) 



{Miip') - B,{p')y 



1 



--{3Io{z) + l2{z))K,{z), 



(78) 



where K ^ 0.29 GeV~^, the instanton density p ^ 1.667 GeV^^, P = \/—p'^, and the 
I's and Ki are modified Bessel functions of the first and second kind, respectively. 

In Ref. 113] a number of calculations were done with both the polynomial infra-red en- 
hancement and the instanton contribution, with various choices of the form to check for 
double-counting. See that references for the results and details on how the calculation is 
done. 

This approach is most promising when coupled with the methods for using SDE solutions 
to get mesonic properties as in Refs. jlHl 1^ HH HH] • In the previous work, in which the 
key idea given by Eq. ()70|) is used, the use of instant form SDE solutions limits the work 
to low momentum transfer, while with the light cone solutions described in this section all 
momentum transfer can be studied. This is an exciting prospect for future research in light 
cone physics. 

6 Deeply Virtual Compton Scattering and Skewed Par- 
ton Distributions 

In recent years there has been a great interest in deeply virtual compton scattering (DVCS) 
arising from the use of this reaction for studying generalized parton distributions |57t 1^ . also 
called off- forward or skewed parton distributions (SPDs). See Ref. [59j for a recent review 
with references. Virtual Compton scattering from a hadron, i.e., the e -|- e' — >^ 7 reaction on 
the hadron, is the reaction 



where q (q') are the momenta of the incoing virtual (outgoing real) photons and P (P'=P-A) 
are the momenta of the initial (final) hadron, H. The DVCS amplitude is 



with T the time-ordered product of the electromagnetic currents. 

DVCS is the high = —q^ limit of VCS, and has some of the physical as well as math- 
ematical aspects of deep inelastic scalleting (DIS). A striking aspect of inelastic scattering 
with a hadronic target H 



j*{q) + H{P) 



7(g) + HiP') 



(79) 




(80) 



e + H 



e' + hadrons 



(81) 
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H 

Figure 12: Handbag diagram for DIS cross section 



is that the differential cross section factorizes 



dE'dn g4 

where L^f^^ is the leptonic tensor 



E'EL'^^'W^,, (82) 



Lf = lTr[(r + me)7^(^ + me)7l, (83) 

with {k^'.k^') the (incoming, outgoing) electron four-momenta, = e^/47rand W^^i, is the 
hadronic tensor. For DIS the "handbag" diagram shown in Fig. dominates. In the 
diagram H is the hadron target and the sum is taken over intermediate states X, expressing 
the hadronic tensor as 

^''^^'^^ = ^|rfVe'^-<i/(P,A)|[JL(x),JL(0)]|i/(P,A)> (84) 

in terms of the matrix element of the current commutator n the hadronic state with spin 
index A. For more than two decades DIS has been used to learn about parton distributions. 
(See, e.g., Ref|Ml) 

For DVCS the amplitude (not the cross section as in DIS) is dominated by the handbag 
diagrams, shown in Fig. ^| In the diagram the target is shown as \E', the B-S amplitude of 
the hadron, which gives the vertex of H, a quark, and the residual system. After a review of 
the kinematics we give a brief review of DVCS on a nucleon and review recent developements 
on skewed parton distributions of the pion, where newly developed B-S amplitudes can be 
used to compare with experiment. 
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Figure 13: Handbag diagrams for DVCS amplitude, a) s-chanel; b) u-channel. 



6.1 Kinematics for DVCS. Skewness 



With the notation for a hght cone four vector = {V^, V , ^l), choosing the frame with 
Pi = 0, and defining the four momentum transfer A = P — P', the kinematics of VCS with 
a hadron target of mass M is 



P' = 



{I -OP' 



(85) 



A = P-P'^ 



with 



(86) 



(87) 



the skewedness parameter describing the asymmetry in plus momentum. The squared mo- 
mentum transfer is 



t = = 2P • A = 



+ Ag 
1-e 



(88) 



In a frame where the incident spacelike photon has g+ = 

<1 = 



^ (q^ + A^)^ , jM^ + Aj ■ 
^' 7" TVTTT' ^-L 



ep+ ■ {i-i)p+ 

(q± + Ax)2 



-,q± + Aj 



(89) 
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In deeply virtual Compton scattering (DVCS) where = —q^ is large compared to the 
mass M and — t, one obtains 



- (90) 

i.e. ^ plays the role of the Bjorken variable in DVCS. For a fixed value of — t, the allowed 
range of ^ is given by 



0<«<fejlf,l.i^-l). (91) 



2M2VV 

6.2 DVCS on the Nucleon 

There has been a great deal of interest in elastic and inelastic form factors of the nucleon for 
many years for tests of quark models and quark distributions. The interest in DVCS was the 
realization that new information on quark structure of the nucleon can be obtained, and it 
can be directly related to the elastic structure functions. First recall the elastic form factors 
of the nucleon, -Fi, F2, 

<n{P')\J';r\n{P)> = N{nb,F,{q')+ta,^£^F,iq')]NiP), (92) 
with N{P) the nucleon spinor, the axial form factors 

<p{P')\Al + tAl\n{P)> = N,{P')[^,^,GA{q') + q,l,HA{q')]NrXP), (93) 
where = N''y^'j5^ is the axial vector current, and the pseudoscalar form factor 

< n{P')\^\n{P) > = N{P')Gp{q')N{P). (94) 
The VCS amplitude with the handbag diagrams. Fig. [121 have the form 

T^- = ^ J 0-Tr[^-uS{k + q)^"'u + ^'^uS{k-q)^''u]J^{k), (95) 
where S is the quark propagator and 

J^{k) = [ e'^^'d^x < N{p')\il){Q)^{x)\N{P) > . (96) 



When the DVCS amplitude is expressed in light cone variables, the quark matrix element 
is expressed in terms of skewed quark distributions which in the notation of Ref . jSTj are to 
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twist two 



2-n 



dz - ~ - 

—e^^-^<N{p')\mi'l^i^{^-)\N{P)> ^ H{x,A',ON{P'h'l,N{P) 

2 .x^..„n75A- 



+Eix,A\ONiP')^iP) (97) 



These skewed distributions satisfy the sum rules 



J dxH{x,A'^,^) = Fi{A^ 
dxE{x,A\^) = F2{A^ 



-1 
1 

dxH{x,A'^,^) = GAiA^) (98) 



-1 
1 



dxE{x,A',0 = Gp{A'). 

In recent years there have been studies of higher twist effects for DVCS. For a recent 
review of DVCS with a nucleon target and references to recent work see Ref. [HS], and a 
hght cone wave function representation of DVCS is discussed in Ref. [^2] ■ There is a great 
interest in new experiments to measure skewed quark distributions. These were discussed at 
the international workshop Light cone 2002 4-. 

We now turn to DVCS and skewed quark distribution on the pion. 

6.3 DVCS on the Pion 

The skewed quark distributions of the pion,jF^(,^, x, t), are defined by the integral 

^e^^^"^-/2(P'|^(0)7+V^(^)|P)|.+=z,=o = J'Atx,t){P + Py, (99) 

where z = {z^, z~ , z_|_) in a light front representation. Note that the path-ordered exponential 
of the gauge field, Pexp[i J z^^A^j], required by gauge invariancefSH] in Eq. does not 
appear in the light front gauge A'^ = 0. Recalling the definition of the pion form factor, 

J+(0) = (P'|7/;(0)7+^(0)|P) = FAt){P + P') + , (100) 

as one can see from Eqs. (jlOO|) and (j99p . the JF^ involves one less integration than the form 
factor due to nonlocality of the current matrix element. The SQDs display characteristics 
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of the ordinary (forward) quark distribution in the hmit of ^ ^ and t ^ 0, on the other 
hand, the first moment of the SQDs is related to the form factor by the following sum 
rules mi EH]: 

/ dxJ^^{i,x,t) = F^{t), (101) 
Jo 

where jF^(^,x, t) = euJ^^{^,x,t) — erfjF^(^,x, t) and we assume isospin symmetry(mu = mj) 
so that jF^(^,x,t) = jF^(^,x, t). Note that Eq. p01|) is independent of ^, which provides 
important constraints on any model calculation of the SQDs. 

The diagrams for the SQDs are shown in Fig. El In the range of x, < a; < 1, for 
X > ^ the BS amplitude is the standard one that determines the pion form fator, as shown 
in Fig. IT^ b). while for x < ^ an analytic continuation, called the nonvalence BS amplitude 
is needed, as shown in Fig. IT^ c). 




Figure 14: Diagrams for SQDs in different kinematic regions for the case ^ > 0: The 
covariant diagram (a) corresponds to the sum of the LF valence diagram (b) defined in 
^ < a; < 1 region and the nonvalence diagram (c) defined in < a; < ^ region. The large 
white and black blobs at the meson-quark vertices in (b) and (c) represent the ordinary LF 
wave function and the nonvalence wave function vertices, respectively. The small white blob 
at the quark-gauge boson vertex indicates the nonlocality of the vertex. 



The B-S amplitudes (see Sec. 4) needed for the calculation of the SQDs are solutions to 

{Ml - Ml)x{x.,^i±) = J[dy][dH^]lC{x,,Kr,yjAj±)x{y,M), (102) 

where /C is the B-S kernel, M| = MV(1 - 0,^o = (^^ + ki)/(l - 2;) - (m^ + ki)/(^ - x), 
and xi^iy^i±) is B-S amplitude. Defining the solution to Eq. ()l()2|l when a; > ^, which 
is the familiar B-S amplitude, as the valence amplitude, the valence and nonvalence B-S 
solutions are 

X""" = X (x < (103) 

(104) 
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In the light front quark modeljnS] the B-S amphtudes in the two regions are related by 

{Ml-Ml)x''''ixi,k,^) = J[dy][d'h]}Cix,,hr,y„l,^)x''''\yj,\,±). (105) 

The contribution from the process shown in Fig lT^ b) is given by the x > ^ B-S amplitude, 
with the valence part of J-'.„{^,x,t) given by 

Nc dx 6{x- k+/P+)AP+ 



Nc f 



IGvr^ xil — x)x'{l — x' 



I rf2k^x™'(a:,k^)(kx-kl + m2)x'™'(x',k'^), (106) 

with the internal momenta of the (struck) quark for the final state given by 

x' = k^ = k^ + i^A^. (107) 

The contribution from the nonvalence part of jF^(.^,x,t), the process shown in Fig lT^ c). is 
given by the x > ^ B-S amplitude, with the form 

^-(e,x,t) ^ 5(x-fc+/P+)4P 



(P + P')+ Jo 167r3 a;(l - x)x'{l - x') 
J d'k^xVix,k'[)H{x,k^)x''\x',k'^) , (108) 



with H{x, kx) a function of kinematic variables given in Refjl6j. x^ is the quark-gauge boson 
vertex shown as the small white blob in the upper vertices of Fig El and k'^ = k± + x" A± 

and X™ 

The contribution from this nonvalence region x < ,^ is obtained by making approximate 
use of the constraint 



FAt) = C dxrA\i,x,t)+ f dxJ'Tii.x.t), 
Je Jo 



(109) 



so the B-S amplitude x^°'\x,'k±) is not used. The valence B-S amplitude is modeled by a 
light cone harmonic oscillator: 

X^'^Ha:,k^) = U— ^ exp(-kV2/32) . (110) 



Substituting Eq. ()llUp in Eq. p06|) one obtains the valence contribution to the SQDs. 

For ^ = the pion form factor is given by the valence process. In Fig^J the valence 
contribution to P^ is shown. For values of ^ > 0.3 the nonvalence contribution is seen to 
become significant. 
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Figure 15: The valence contribution to the pion EM form factor with different skewedness 
parameters compared with the experimental data 



In Figs. IT^ and IT71 are shown the SQDs, t), of the pion for fixed momentum 

transfer -t = 0.2 GeV^ (0 < ^ < 0.92) and -t = 1.0 GeV^ (0 < ^ < 0.98) but with differ- 
ent skewedness parameters ^, respectively. The solid and cross (x) lines in the nonvalence 
contributions are the exact solutions obtained from using the B-S equation to obtain ^™ for 
the model of Eq. ljllOj) . and an approximation based on Eq. p09j) . The SQDs at ^ = as 
shown in Figs- ITBT a) and lTTT a) correspond to the ordinary quark distributions with vanishing 
nonvalence contributions. The frame-independence of the model calculation is ensured by 
the area under the solid lines(valence + nonvalence) being equal to the pion form factor at 
given —t. As one can see from Figs. ITHT b-c) and lTTT b-c). while the nonvalence contributions 
are small for small ^ = 0.3, they are large for large skewdness parameter ^ = 0.9. It is 
interesting to note that the instantaneous contribuitons (dotted lines in each figure) become 
more pronounced as ^ — >■ ^max for each —t. While the instantaneous part of the nonvalence 
contribution vanishes as x — ^~ = lim^^Q{^ — e) as shown in Figs. El and El the net result 
of J-'™{^,x,t) including the on-mass shell propagating part does not vanish as x = ^ and 
consequently causes a discontinuity to the zero value of ^™'(^,^,t). However, such discon- 
tinuity at a; = ,^ is just an artifact due to the difference in the x —>■ C, behavior between the 
gauge boson vertex {x^ in Ref.|16j) and the hadronic vertex (x^"' (^')kj_)) in Eq. p06p of 
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our approximate model calculation. See the following subsection for a discussion of and 
the continuity problem. We have indeed confirmed that the discontinuity at x = ^ does not 
occur in the limit of a point hadron vertex as observed in the QED calculation j,62j. Thus, 
for a full analysis of DVCS satisfying factorization theorems, it would be necessary to solve 
the bound-state B-S equation similar to Eq. ()lU2p for the gauge boson (x^) as well as for 
the hadron (x™')- 

The main approximation in this calculation is the treatment of relevant operator /C(x, k^; y, l±) 
in Eq. ()105|) connecting the one-body to three-body sector by taking a constant G-,, for the 
quantity G-,, = f[dy][(f\±]]C{x,'k±;y,\±)x^°'\y,i±), which in general depends on x and kj_. 
The reliability of this approximation was checked by examining the frame-independence of 
the numerical results and using the sum rule given by Eq. (|1(J9|) . This method seems useful for 
the present study of the relation between SQDs and the form factor in the nonperturbative 
regions. 



6.3.1 Continuity of Pion DVCS Generalized Parton Distributions 

The problem with continuity, discussed in the previous section, is the treatment of the gauge 
boson vertex, iii the nonvalence region, ^ > x; i.e., the white blob in Fig lT^ c). In RefjTB] 
the form used was 

x\-"^k:) = (111) 

with x" = x/x, k'[ = k_L+x"A± and M"^ = {k"^ +w?)/x" + {J^'^ +m?)/ {I - x"). In [T7I an 
effective light cone wave function was used for the gauge boson wave function: 



ViV,/3x"(l-x") 



This model solves the problem of continuity, as shown in Fig. UHl 
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Figure 16: Skewed quark distributions of the pion at — t = 0.2 GeV^ with ^ = in (a) 
,0.3 in (b), and 0.9 in (c), respectively. The sohd [cross (x)] hne in nonvalence contribution 
represents the full result of using true [average] value and the dotted line represents the 
instantaneous part of the nonvalence contribution. 
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Figure 17: Skewed quark distributions of the pion at — t = 1 GeV^ with different ^ = in 
(a), 0.3 in (b), and 0.9 in (c), respectively. The same notation is used as in Fig. [TBI 
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Figure 18: Generalized parton distributions of the pion with modified gauge boson wave 
function at —t = 0.5 GeV^ with ^ = 0.3 in (a) and 0.6 in (b). 

7 Rare B Ki+i' Decays 

There is currently a great interest in the study of B decays for study of the CKM mixing 
matrix as a test of electroweak models. The exclusive B Ki~^i~ decay, a rare decay in 
that it is forbidden for tree-level processes, is particularly interesting, but the theoretical 
calculation requires hadronic form factors involving nonperturbative physics. There have 
been many theoretical studies in recent years[nni EIlinHlinnilZnilZIllZaiZSllZllZSllZni- 
Since this decay of B-mesons involves a large energy transfer, it is a natural for a light cone 
treatment, which has many advantages compared to instant form relativistic quark model 

treatments inni EH ESI EH] • 

The effective Hamiltonian for the b si~^i~ decay process in obtained after integrating 
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out the heavy top quark and the bosons [77j : 

AG 



= -|\4,\/,:5^a(/x)o.(^), (113) 

* i 

where is the Fermi constant, Vij are the CKM matrix elements and Ci{fi) are the Wilson 
coefficients. The important operators Oi{fi) to the rare b —>■ s£~^£~ decay are 



Oi = 




02 = 




07 = 


-r^mk{s^a^^PRb^)F'"' , 

lOTT^ 


09 = 


t|^(s.7'^Pl&.)(^7/), 


Olo = 


^^{sarPLhMl.ld), (114) 


1 -F 75)72 is the chiral projection operator and F'^" is the electromagnetic 



where Pl{r) : 

interaction field strength tensor. The Wilson coefficients Ci{m}j) determined by the renormal- 
ization group equations(RGE) from the perturbative value Cj(M^y) are given in the hterature 
(see, for example [751 l79j). After integrating out the terms involving cc- loops and neglecting 
the strange-quark mass, the resulting effective Hamiltonian corresponding to Eq. (|113|) is 



v/2 167r2 



+ Cf{m^Ys^^PMin + Cr^{m^Ysi^PLUi^id 



(115) 



where the effective Wilson coefficient Cf{s=q^lml) is given in Refs [ZHl EHl IHOl IHD ISS ES! , 
and has the form 

cf{s) ^ c'f(s) + rLD(s), 

= C9(^l + ^a;(S)^ +rsD(5)+lLD(S), (116) 

where the functions 1sd(s) is the one-loop matrix element of O9, 1ld(s) describes the long 
distance contributions due to the charmonium vector resonances via B —>■ K{J/ip, ip',- ■ ■) ^ 
Ki~^i~ , and a;(s) represents the one-gluon correction to the matrix element of O9. 

The parameters used here, given by Refs. |H2l ES], are mt = 175 GeV, = 4.8 GeV, 
rric = 1.4 GeV, asiMw) = 0.12, a,(m6) = 0.22, Ci = -0.26, C2 = 1.11, C3 = 0.01, 
C4 = -0.03, C5 = 0.008, Cq = -0.03, C7 = -0.32, Cg = 4.26, and Cio = -4.62. 

The long-distance contribution to B ^ K decay is contained in the meson matrix ele- 
ments of the bilinear quark currents appearing in Ties- The matrix elements of the hadronic 
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currents for B K transition can be parametrized in terms of hadronic form factors as 
follows: 



= {K\sYPLb\B) 
1 



(117) 



and 



= {K\sia^"'q^Pnb\B) 



;ii8) 



where P = Pb + Pr and q = Pb — Pr is the four-momentum transfer to the lepton pair and 
Amf <q^ < {Mb - Mr^. 

The transition amplitude for the B Kt^t' decay can be written as 



M 



{Kti~\n,^\B) 



4:Gf Oi 
72 47r 



;ii9) 



where a = e'^/An is the fine structure constant. The differential decay rate for B — > Ki'^i 
with nonzero lepton mass(m^ 7^ 0) is given by 



dr 

ds 



X 



(120) 



where 



\CfF^ - 



2C7 



rp |2 I |2|7Ti |2 



1 + Vr 

Cion(l-r)VoP-0|i^+P] 



= (s — 1 — f) — 4f, 
with s = q^/Ml, rhi = mj/M'^, f = M]^/Ml, and 



(121) 



2\ I q F ^„2^ 
Ml -M2 



(122) 



Note fEqs. (|12Ull2ip ) that only -F+(g^) and Fxi^q^) are necessary for the massless(mf = 0) 
rare exclusive semileptonic h — > s£+£~, with Fo+{q^) not contributing. 
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The differential branching ratio, dBR{B — > Kl^C, )/ds is obtained from 
dBR{B^Ki+£-) dr{B^Ki+i~) 1 , 



ds 



(123) 



tot 



G \ 

with Ftot = g|^3^ v^ch\ the total width. Another interesting observable, the longitudinal 
lepton polarization asymmetry(LPA) is defined as 



Pl{s) 



dTh=-i/ds — dTh=i/ds 
dTh=_i/ds + dTh=i/ds 



(124) 



where h = +1(— 1) denotes right (left) handed i in the final state. From Eq. p20|) . one 
obtains for B Ki+i- 



Pr.(s) 



2^^1-4^ 















(125) 



The calculation of i? — * K£~^i~ is similar to the calculation of !Ftt{^, x, t) described in the 
previous section. An important difference is that in calculatin F^{q'^) and Fxiq^) the initial 
and final B-S amplitudes are different, so that in contrast to Eqs. ()106|l and ()108|1 with x 
the initial and final state, one needs a xi ^^id X2 for the initial b-meson final s-meson. 

First, if one works in g"*" = frame, with = —q]_ < Q the weak form factors needed are 
of the form 



7^ ■ 



dx / d'^k±Hi{x,k±,k'j_)(j)2{x,k'j_)(f)i{x,kj 



and 



Friq 



dx / d'^k±H2{x, k±, k'j_)(f)2{x, k'j_)(l)i{x, k±) 



(126) 



(127) 



with the known kinematic functions Hi and H2 given in Ref |15j . The model light front 
Gaussian wave functions used are 



{x, k±) 



1/2 



exp{-ky2p^ 



;i28) 



normalized by J d^k\(j){x, k±)\'^ = 1 where k"^ = kj_ + k1 and kz is defined as 



1, 



2 2 
— m| 



'^-2'*'"+ 2M„ 



(129) 



The calculations in Ref|15| are done in the purely longitudinal frame with g"^ > and 
Pix = P21. = 0, and thus the momentum transfer square (f' = q^q^ > is timelike. For the 
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valence process, corresponding to Fig lT^ b). the calculation is carried out using the model 
1-f wave functions k±). For the nonvalence process, corresponding to Fig lT^ c). the final 
state B-S amplitude must be replaced by the function 

GBKix,k±) = [ [ (fl±H3{{x,k±;y,l±)IC{x,k±;y,l±), (130) 

where is a known function of the kinematic variables. As discussed in the treatment 
of skewed quark distributions in the previous section, in Ref [T3] Gbk{.x, k±) is taken as 
a constant, and the reliablity of this approximation is checked by examining the frame- 
independence of the numerical results. 

The results with this model for the differential branching ratios, with = and 
mr=l.777 GeV, for B Ki'^i'li = e,/i) are shown in Fig. fTW a) and for B Kt^t~ 
in Fig. ITW b). respectively. The thick(thin) solid line represents the result with(without) 
the LD contribution(FLD(s)) to Cf given by Eq. pi6|) . One can see that the pole con- 
tributions clearly overwhelm the branching ratio near J/ip{lS) and ip'{2S) peaks, however, 
suitable invariant mass cuts can separate the LD contribution from SD one away from 
these peaks. This divides the spectrum into two distinct regions jHH EH] : (i) low-dilepton 
mass, imj < < Mj^^ — 6, and (ii) high-dilepton mass, M|, + S < q"^ < Q'^ax^ where 
6 is to be matched to an experimental cut. The numerical results for the non-resonant 
branching ratios(assuming \Vtb\ ^ 1) are 4.96 x 10~'''|Vts/V'cfeP for B Ki^i~ {£ = e,/x) 
and 1.27 x 10-^|\4,/KbP for B Kt+t', respectively. While the data by CLEO Col- 
laboration jHEl reported the branching ratio Bt[B — >■ Ke^e~) < 1.7 x 10~^, the data 
by Belle Collaboration(K. Abe et al.) ^ reported Br{B iTe+e") < 1.2 x 10"^ and 
By{B Kfi+^') = {0.99toiltoil) X 10"*^, respectively. The comparison of the results of 
this model with other models is given in detail in Ref J3] 

The longitudinal lepton polarization asymmetry(LPA) in this model for B Kn'^fi~ and 
B Kt^t~ as a function of s, respectively, in Figs. EUT a) and (b). The results with the LD 
contributions are shown be the thick solid line, and without by the thin solid line. Further 
discussion of these results and comparison with other calculations is found in Ref J3] . 

In conclusion, the rare exclusive semileptonic B Ki'^i" {i = e, fi and r) decays are 
found to be useful for studies of nonperturbative form factors, and the light cone approach 
is most valuable for the theoretical models. 
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Figure 19: The branching ratios for B Ki^i'la) and B Kr^r^ih) transitions. The 
thick(thin) sohd hne represents the result with(without) LD contribution(yi£)) to Cg^ in 
Eq. (jl2H) . The dotted hne in (b) represents the FqIq"^) contribution to the total branching 
ratio of r decay. 
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Figure 20: The longitudinal lepton polarization asymmetry Pl{s) for B 
B KT^T~{b) transitions. The same line code is used as in Fig. IT^ 
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8 Meson Light Cone Wave Functions 



The structure of mesons is of great interest for the study of QCD. With a simple two- 
body configuration qq being the lowest Fock state, for certain reactions the pion can be 
a relatively simple system for theoretical investigations. As discussed in Sec. 4, for high 
momentum transfer processes the PQCD treatment is adequate. At low and medium mo- 
mentum transfers nonperturbative effects dominate. In the present section the QCD sum 
rule method for determining the wave function of the tt and p mesons is reviewed. The pion 
wave function is of particular interest as there are now experimental studies of the light cone 
wave function, which is discussed in this section. 

In Sec. 4 the Bethe-Salpeter amplitude for a two-body bound state, '^{zi,Z2) = 
< 0|T[^(2;i)^'(2;2)] |M > in instant form coordinate space, was discussed, with the light cone 
B-S amplitude, '^{xi,X2,k±), given as the solution to the 1-c B-S equation, Eq.( HT]). In 
order to obtain the correct physical B-S amplitude, also called the light cone wave function, 
one needs a kernel {K{xi,X2, k±,yi,y2,l±) which gives a satisfactory representation of the 
field theory being considered. In the applications discussed in previous sections models were 
used. In the next subsection we discuss the use of QCD sum rules to estimate the pion and 
rho meson wave functions. Theoretical and experimental work on the direct measurement 
of the qq component of the pion wave function is discussed in the following subsection. 

8.1 QCD Sum rules for the Pion and Rho Meson Wave Function 

We first discuss the pion wave function is some detail, and then give the results for the rho. 
The light cone wave function of a tt"*" meson can be obtained from the gauge-invariant matrix 
element 

^f.{z-p,z^) = <0\d{0)j^Ve^ody'^-j,u{z)\7r^ip)> , (131) 

where (u,d) are the up and down quark fields and = A^X"-/2 is the QCD gluonic color 
field. We shall work in the fixed point gauge, with the gauge condition x^A^ = 0, so that 
the operator e-^o "^^^^m = i. xhe method of QCD sum rules[87j utilizes an operator product 
expansion, which leads in a natural way to a twist expansion of the wave function. We briefly 
review QCD sum rule procedure 

8.1.1 Brief Review of the QCD Sum Rule Method 

The method of QCD sum rules makes use of a correlator defined in terms of a composite 
field operators rja, rjb 

Uabip) = I j (fx e^'-" < I T[r/,(x)r/,(0)] | > . (132) 

For mesons the ri{x) can be taken in the form 

Va{x) = q{x)Taqix) , (133) 
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with Ta chosen so that the operator rja creates states with the quantum numbers of the 
meson under consideration. E.g., for a scalar, pseudoscalar or vector meson Ta = 1, 7^x75 or 
7^, respectively, are possible choices. 

If the correlator satisfies an unsubtracted dispersion relation (in Euclidean space), then 
Uab can be expressed as 

Hab p = - / ds 

where M is the mass of the stable meson. Nab is a constant, and So is the parameter for the 
beginning of the continuum contribution to the dispersion relation. The QCD sum rule is 
obtaind by equating the Borel transformation, i3_p2^Af|) of the phenomenological correlator, 
Eq. p34|) . to the QCD evaluation of liab^ 

n,,(Ms) = iV,,e-*^'/^^ + i / dse^^'^^Bimliabis) 

^ J So 

= ^T^'iMs), (135) 

with Mb called the Borel mass. The Borel transformation reduces the importance of the 
continuum contribution, and results in a rapidly diminishing size of the QCD diagrams with 
increasing dimension. The QCD side of Eq ]135l is of the form 

n?r(MB) = K^'^'iM,) + ^ , (136) 

where the c„ are constants obtained by evaluating the QCD diagrams and the < 0„ > are 
vacuum condensates of dimension n. In order of increasing dimension the most important 
three condensates are 

< O3 > = < qq > quark condensate 

< O4 > = < {C^")'^ > gluon condensate 

< O5 > = < qa^^G^^q > mixed condensate . (137) 

The quark condensate is rather accurately known, with — (27r)^ < qq >~ 0.55 GeV'^. A 
typical value for the gluon condensate is < g^G^ >~ .47 GeV^, and for the mixed condensate 
< gsqfy^^'G^uq >= -ml < qq >, with ml ^ 0.8 GeV^. 

8.1.2 QCD Sum Rules For the tt Light Cone Wave Function 

A QCD sum rule analysis of \E'^(x-p, x^) to extract an expansion of the pion 1-c wave function 
with nonperturbative QCD effects was carried out in Ref. f^. See Refs. |5H1 IH^ for reviews. 
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Figure 21: Diagrams for QCD calculation of 0^ up to D=6. (a) is the PQCD process, (b) is 
the D=4 gluon condensate process, (c), (d) are the processes with two quark condensates 



A study of the 1-c wave function vs. light front quark models has also been done using QCD 
sum rule techniques jHO] 

Using the fixed point gauge and expanding u(z) about the origin, one obtains for 
(Eq.dlSID) 

<0\d{0)^^^,u{z)\7r^{p)> = S„^-^ < 0|d(0)7^75(^.I^")l7r(p) > 

= pu<pTv{z ■ p,z^), (138) 

. with D'^ = d'^ — igA'^ and proportional to the 1-c wave function. 

This expression is evaluated using a version of QCD sum rules. Using t]tj+^{x) = 
d{x)'jfj,'-f5u{x) for the pion composite field operator, the n*^ term in Eq. p38j) is given by 
the vacuum matrix elements 

^"(P') = -p^^^ j dxe'P-^ <0\d{x) ^^.(z-nruixMO) M0)\0> . (139) 

By evaluating the diagrams up to dimension 6, illustrated in Fig.|^ Chernyak and Zhitnitsky 
derived the expansion for the Borel transform of J„(p^), 

^"^M^) = 7r^4-^ + -«^f^4S^ + ^<0|(2vr)^gg|0>^V---a40) 



The light cone wave function is obtained from (l>n{z ■ p, z^) of Eq. ()138|l by letting 7^ 
7+ = (7*^ + 7^)/a/2 and z'^ = z± = 0. The light cone wave function is then obtianed from 
\E'(2;~p+, 0) in the gauge with = : 



rl7~ 

/ ^e'"""^'+^(2-p+,0) . (141) 
Joo 27r 
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The CZ wave function |24j. with the normahzation J dx(j){x) = with the pion decay 
constant = 133 MeV, was obtained using the values of the condensates given above. 
Within the errors in the method an approximate form for the CZ result often used is 

0^^(Ob=.5GeV = 5x{l-x){l-2xf , (142) 

with units of U/ ^3 /2) This can be compared to the asymptotic light cone wave function j27j. 
also with units of f^,/ a/(3/2) 

<Pf{x)\Q2_,^ = x{l-x). (143) 

8.1.3 QCD Sum Rules For The Light Cone Wave Function of the p 

Treatment of the rho meson, with the composite field operator t]^ is quite different for the 
longitudinal {tj^^q) and the transverse {r]^^_^i) currents, for and p±, respectively. 

For pl the QCD sum rule treatment of Ref sM, is the same as for the tt, with the 
quantities Inip^) given by the expression in Eq. ()139|) with the operators 75 removed. The 
resulting expansion of the 1-c wave function, using the analog of the diagrams in Fig. |^ in 
units of 3fp is approximately 

02(OlQ=.5GeV = x(l -x)(0.7 + 1.5(1 - 2x)2) , (144) 

For the quantities corresponding to the Inip^) given by the expression in Eq. (jl39|) for 
the TT are obtained by replacing Zy'-)^'~^^ by a^^Zy in Eq. (|139|) (with a sum over p, of course). 
The result of the analysis [24j for the transverse rho 1-c wave function, in units of 3fp is 
approximately 

02(OlQ=.5GeV = 5(x(l-x))% (145) 

In Fig. 1221 the CZ pion 1-c pion wave function is compared to the asymptotic 1-c pion 
wavefunction, in units of /^/ a/3/2, for the pion and the transverse rho 1-c CZ wavefunction 

9 Factorization and Experimental Study of Light Cone 
Wave Functions 

An essential aspect of testing hadronic and nuclear structure by high energy reactions is 
that amplitudes for such reactions factorize. Consider Fig. |2S1 with a schematic high energy 
reaction on a target with wave function 0j„ leading to a final state with wave function 
(pout- The factorization theorem states that the reaction amplitude satisfies the schematic 
relationship 

M{in,p out,p') = dx ddy^in{x,Pi)T{x,y,Q'^)(l)out{y,Pf) , (146) 

Jo Jo 
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Figure 22: Comparison of the CZ with the asymptotic tt hght cone wave function, in units 
of /^/-\/2, and the transverse p CZ 1-c wave function. 




Figure 23: High-energy reaction of a projectile on a target (in) leading to a final state (out). 
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Figure 24: Coherent diffractive dissociation of a pion projectile into two jets with a nuclear 
(X) target 



with 0j„, (pout the initial and final target wave functions, and T represents the parton hard 
scattering. Factorization is important for the extraction of information for many reactions. 
For a discussion with references see, e.g., Ref. ||89j. The coherent production of dijets with 
a high energy pion projectile, which we discuss next, is of particular interest because of a 
recent experimental measurement. 

9.0.4 Coherent Diffractive Pion Production of Dijets 

The concept of diffractive dissociation is very old in nuclear/particle phvsics^^. and was 
extended to QCD in Ref. [HH]- The essential idea for pion production of dijets is that for a 
pion projectile on a nuclear target at high energy, the short-range color zero qq configuration 
will traverse the nucleus with little interaction, while the other components will be absorbed. 
Therefore the qq component is filtered out and leaves the nucleus to form two jets. 

The (tt, dijet) amplitude [Sll, for high energies with the forward qq scattering amplitude 
f{b^) ~ sa{b'^), where cr(6^) is the total cross section as a function of impact parameter b, 
with this coherent diffractive dissociation production scenerio is 



where b is the impact parameter, a{b) is the total forward cross section of the valence 
diquark qq system, and the scattering of the final jets has been neglected, so that their 
wave function is a plane wave with the relative transverse momentum variable k±. This is 
a special case of the factoriation theorem, Eq. ()14(i|l . Recently, there has been a theoretical 
study of photoproduction of high energy dijets with nuclear targets[95] which within the 



M 




(147) 
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approximations leading to Eq. (jl47|) is the reaction 



7 + TV (qq) + TV , 

so that the pion valence wave function can also be studied with this reaction. 



(148) 



9.0.5 Experimental Study of the Pion Light Cone qq Component 

A direct measurement of the light cone meson wave function can be carried out through 
a high-energy coherent diffractive dissociation of the meson. This was recently done at 
Fermilab experiment E791|ni]. The experiment used a 500 GeV vr^ beam with a platinum 
target. The diffractive dijets were identified as having their yield depend on the transverse 
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Figure 25: x distribution of dijets in theE791 TV(7r, dijet)TV reaction. The dashed line is 
the asymptotic, the dotted line the CZ and thesolid line is a combination of the CZ and 
asymptotic 1-c wave functions. 



momentum transfered to the nucleus, by 



jet 



-<b'^>ql 



(149) 

where < 6^ >= R%-/3 is the effective impact parameter squared. The x value of the jet (see 
Fig. is obtained from 



X. 



Pjetl 



exp 



Pjetl + Pjet2 
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(150) 



The results of the experiment are shown in Fig. 123 The fit to the data, the sohd fine in 
Fig 1221 is made with a squared wave function that is a hnear combination of the CZ and AS 
wave functions 



<P\x) = aAsK {xY ^ acz<K\^r ■ (151) 
Using the relationship for the virtuality of the dijets 

- (152) 

for A;_L > 1.5 GeV >~ 10 GeV^, and the ^as gives a good fit to the data, while in the 
lower 1.25 < 1.5 bin the CZ component makes a sizable contribution. 



10 Conclusions 

The methods of light cone field theory and light front quantum mechanics play a vital role 
in the extraction of information on the structure of hadrons and nuclei. It has been applied 
in recent years to hadronic form factors and a variety of reactions, with a great deal of work 
on the transition from nonperturbative to perturbative QCD treatments of hadronic form 
factors and deeply virtual reactions with photon projectiles. The use of light cone wave 
functions is important for extracting information from the decay of heavy-quark hadrons. 
Quite recently coherent diffractive dissociation and related reactions have begun to provide 
direct measurements of light cone wave functions. This is an exciting field for theory and 
experiment, and will continue to provide important informationabout the nature of Quantum 
Chromodynamics. 
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